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The magnetic properties of a variety of materials with promise for technological applications 
have been described by models in which fermions are coupled to local moment spins. Monte Carlo 
studies of such models usually ignore electron-electron interactions, even though the energy scale 
corresponding to these interactions may be comparable to or larger than other relevant energy scales. 
In this work we add on-site interactions between fermions to the double exchange model which we 
study with a Monte Carlo scheme in which temporal fluctuations of local moment spins are fully 
accounted for and electron-electron interactions are treated at a mean field level. We show that 
when the number of fermions is considerably less than the number of local moments even moderate 
interactions can lead to significant enhancement of ferromagnetism and the Curie temperature. 
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The magnetic properties of a variety of materials of 
fundamental interest and technological promise, such as 
collosal magnetoresistance (CMR) manganites [II], rare 
earth hexaborides [2[, and diluted, magnetic semicon- 
ductors (DMS) [3j have been described using double ex- 
change (DE) type models in which fermions are coupled 
to local moment spins. In the limit that the number of 
fermions is considerably less than the number of local 
moments, such models generally display ferromagnetism 
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In manganites it is well established that electron- 
electron interactions are at least as large as the Hund 
coupling, and may be the largest energy scale in the 
problem [l|, l5|— |tJ . The importance of interactions has 
also been stressed for DMS [8l-ll0l] and hexaborides [11] . 
The combination of disorder and electron-electron inter- 
actions may also play an important role in nanoscale elec- 
tronic phase separation |12{ , which has been argued to be 
important for CMR in manganites 13]. It is hence im- 
portant to develop accurate techniques that can account 
for the effects of electron-electron interactions in DE type 
models and to study their effects on magnetic properties 
of these models. 

Here we introduce a method to study the effects of 
electron-electron interactions in DE type models when 
the energy scale for interactions is in the experimentally 
relevant regime of no more than a few times the Hund 
coupling. We use this method to show that even mod- 
erate interactions can enhance ferromagnetism and the 
Curie temperature T c significantly when the number of 
fermions is considerably less than the number of local 
moments. 

We combine a Hartree-Fock treatment of electron- 
electron interactions with a Monte Carlo scheme for 
fermions coupled to classical local moment spins. 
Through the use of exact diagonalizations on small sys- 
tems we have determined regions of parameter space 



where this technique should be most accurate. We calcu- 
late the magnetization as a function of temperature for a 
variety of interaction strengths and determine the effects 
of interactions on the Curie temperature T c . Electron- 
electron interactions can lead to ferromagnetism in the 
absence of Hund coupling 14], so the enhancement in 
magnetization that we find can be understood as this 
tendency reinforcing the ferromagnetism that arises from 
the Hund coupling. 

The general Hamiltonian we consider is of the form 
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where i and j are site indices, iy is the hopping integral, 
Jij is the Hund coupling, U is the on-site Hubbard repul- 
sion, c ia creates a fermion with spin a on site i, Si is a lo- 
cal moment spin on site i, and Sj = g Yl a 8 ( c 'ia <T( ^P c 3? 
is the fermion spin on site j. For simplicity, we as- 
sume the hopping is to nearest neighbour sites on a cu- 
bic lattice, and that the Hund coupling is purely local: 
Jij = J5ij. We assume that there is a local moment 
on every site in the lattice, and that the fermions have 
spin-i (these assumptions can be easily relaxed). 

There has been some work to study the effects of finite 
U in models of the form Eq. ([1]) using exact diagonaliza- 
tions [15[ , DMRG [16[ and mean field approximations [9( . 
However, these techniques are not appropriate for study- 
ing finite temperature magnetic properties in dimensions 
higher than one taking into account the temporal fluctu- 
ations of spins. In order to determine the finite tempera- 
ture magnetic properties of the Hamiltonian Eq. pj). it is 



necessary to perform Monte Carlo simulations. If there 
are N local moment spins with spin S, then the size of 
the Hilbcrt space scales as (25 + 1) N , hence it is usual to 
approximate the local moment spin as classical, which is 
often reasonable given that in many systems of interest 5 
is larg er than \. Previous such Monte Carlo simulations 
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have restricted their attention to models of the 
form Eq. dl]) with U = with the exception of Ref . [8] . 
In Ref. |8| interactions were included for some parame- 
ter values using a zero temperature variational procedure 
due to convergence issues with Hartree-Fock and did not 
appear to have a strong influence on magnetic proper- 
ties. The model in Ref. [8| also included disorder and 
spin-orbit coupling, and these terms may have influenced 
the convergence of Hartree-Fock calculations. 

We write the classical local moment spins in the form 
Si = (Sf,(f>i), and for a specific arrangement of local 
moment spins, one can diagonalize Eq. (p} with U = 0: 



n({S l }) = J2E m ({S l })at 



(2) 



where {E m ({Si)} are the eigenvalues for a given local 
moment spin configuration. a/ m and a m are the creation 
and annihilation operators for the m th eigenstate of H: 

a m = / J VimtjCjcr, CL m — / J Wima C i<T- 



This allows one to write the fermion free energy as 

27V 

^{ s ^)--^E ln ( 1 + e "' (£m({Sl}) ^ ) )' ( 3 ) 



p 



where /3 = 1/fcsT is the inverse temperature and hence 
the classical partition function for the N classical spins 
takes the form 



N 

n 



dSf 



3 -/3^({S,}) 



(4) 



Casting the partition function in this form makes it clear 
that we may use the Metropolis algorithm to determine 
whether or not to flip a spin, with the change in the 
fermion free energy determining whether a spin flip is 
accepted or rejected. 

Introducing a finite U greatly increases the size of 
the fermion Hilbert space: for N local moments and 
n fermions there are (2N)\/(2N — n)\n\ states, as com- 
pared to 2N non-interacting states. This renders even 
moderate values of N out of reach computationally. In 
order to explore larger values of N the interaction term 
must be treated in an approximate fashion, hence our use 
of Hartree-Fock as the simplest self-consistent approach. 
The Hartree-Fock approximation reduces the size of the 
fermion Hilbert space to 2N and allows for the use of the 
Monte Carlo scheme outlined above, with Hartree-Fock 
energies replacing the E m ({Si}). 



We decompose H mt in the Hamiltonian Eq. ([T]) as 
U Y^ nifmi ~ U^2 [( n ^) 4i^l + ( n *l) 4r c *T 

i 

c lt c 4/ 4i c n - (4i c n) 4fCn. 
- (n* t ) K+) + (4t c y ( c li c «t) -( 5 ) 
Approximating _ff int with Eq. ([3]) we can write 

H ~ H HF -uJ2 [( n n) M - (4t c a) (4i c n) . 

i 

and the single particle states satisfy 

HHF \4>rn) = Cm \4>m) , 

with single particle energies 

i 

c !t C *i) ^imlAmf ~ \4l c *t) IPim^iml » ( 6 ) 

which we calculate by iterating to self-consistency and 
use to determine the approximate fermion free energy 
for each local moment spin configuration {S^}. 

The Hartree-Fock approximation is uncontrolled, 
hence in order to ascertain the regions of parameter 
space in which the hybrid Hartree-Fock Monte Carlo 
scheme should be most accurate, we performed exact di- 
agonalization calculations on 8 site 2x2x2 systems 
with 2, 3, and 4 electrons, and compared the energies 
of the exact ground state, first excited state and sec- 
ond excited state averaged over 25 different local mo- 
ment configurations. In Fig. Q] we show the relative error 
A £ = \E n - E% F \/\E n \ as a function of J ft and U/t , 
where E n is the exact energy and E^ F is the Hartree- 
Fock approximation to the energy. 

Provided U is not too large in comparison to J, the 
Hartree-Fock approximation gives a good account of the 
low- lying energy levels, and we found that for U < 3 J the 
relative error was less than 5% with the error decreasing 
with decreased U/J, hence we restricted our Monte Carlo 
simulations to this region. We additionally checked the 
relative error in the densities between the exact diagonal- 
ization results and the Hartree-Fock approximation and 
found the best agreement in the same region of parameter 
space as for the energy. 

We performed Monte Carlo simulations for cubic sys- 
tems with N — 4 3 = 64 and N = 6 3 = 216 local mo- 
ments, with n — 8 and n — 27 fermions respectively 
corresponding to n/N = 1/8 in both cases. These sys- 
tem sizes are competitive with recent simulations in non- 
interacting syste ms [ 20l| . We used a similar equilibration 
procedure to Ref. [19( , by bringing two replicas of the sys- 
tem into equilibrium, evolving under Metropolis dynam- 
ics, with an additional self-consistent Hartree-Fock loop 



as discussed above. We used the z-test 2l| to determine 
when there was 95% confidence that the two replicas were 



in equilibrium and then collected data for a further 10000 
Monte Carlo sweeps. 
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FIG. 1: Relative error between the exact and the Hartree-Fock evaluation of the ground state, first excited state and second 
excited state energies for an 8 site system, with 2, 3 and 4 electrons, averaged over 25 local moment configurations, as a function 
of U/t and J/t. 



We calculated the magnetization of the local moments tion, we use the Binder cumulant [2c 



M (T) = ^ 
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as a function of temperature, where the angle brackets 
indicate a thermodynamic average. In Fig. [5] we show 
the local moment magnetization when J/t = 5 calcu- 
lated both for N = 64 and N — 216 local moments. 
There is some enhancement of the magnetization with 
increasing U/J in the TV = 64 samples. However for the 
larger N — 216 system the enhancement is much clearer 
and is more meaningful as the N — 216 data should 
have smaller finite size effects and be more reflective of 
the thermodynamic limit. Calculations of the fermion 
magnetization yield similar enhancement with increasing 
interaction strength 22 1. 

Due to the evident finite size effects in the magnetiza- 
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to determine T c as a function of interaction strength, as 
g(N, T) should be independent of system size at T c . Our 
results for J/t = 5, determined using N = 64 and N — 
216 are displayed in Fig. [3] There is a monotonic increase 
in T c with increasing U/J, with an almost 50% increase 
in T c between U/ J = and U/ J = 2. 

By performing Monte Carlo simulations of an inter- 
acting DE model, we have demonstrated that electron- 
electron interactions, which are known, or expected to 
be important in a variety of materials whose magnetic 
properties are described by DE type models, can lead 
to quantitatively important increases of the Curie tem- 
perature^ The Hubbard model can display ferromag- 
hence it is natural that in the presence of 
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netism 

Hund coupling, which independently promotes ferromag- 




0.01 



0.02 



0.04 



0.08 



T/t 



FIG. 2: Local moment magnetization for J/t = 5 for N = 64 
and 216 local moments for U/J = 0.0, 0.5, 1.0, 1.5, and 2.0. 
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FIG. 3: Curie temperature as a function of U/J for J/t — 5. 



netism, the two terms in the Hamiltonian can reinforce 
each other to produce the enhancement observed here. 
As a mean field approximation, Hartree-Fock will almost 
certainly overestimate the tendency to ordering as U/J is 
increased, but the effects observed here should be robust 
to fluctuations. This is because the interactions enhance 
an existing tendency to ferromagnetism rather than im- 
posing order on an otherwise disordered system. Fur- 
ther, our use of Hartree-Fock is confined to values of U/J 
where exact diagonalizations suggest that it will be most 
accurate. 

The scheme we have introduced, and our demonstra- 
tion that it should be accurate for the experimentally 
relevant interaction range U < 3J, should be a further 
step towards the quantitative description of the magnetic 
properties of important DE type materials. It should 
also be possible to extend the scheme introduced here 
to treat both disorder (which has been included in pre- 
vious non-interacting simulations [8|, Il8l420| ) and long- 
range Coulomb interactions, which have been argued as 



being relevant to nanoscale phase separation and CMR 
.12]. 
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